! ***************************************
/Finish
/clear
/filn, E52 
/title, Buckling analysis of short cylinder
/units,si	!ùʵλ
! ********* 趨*********
Di=1000	! Ͳо
t=10	! Ͳ
L=5000  ! Ͳ峤
kmodel=2	! ģͿأ1ȫģͣ21/2ģ
nhoop=20	! Բܷ߶εʷ
*if,kmodel,eq,1,then	! ģ趨Ͳʷ
    naxial=40
*else
    naxial=20
*endif
nt=1       ! ȷ߶εʷ
p=1	   ! λѹ
EXX1=2e5	! ϵĵģ
NUXY1=0.3	! ϵĲɱ
!****************ǰ***************************
/PREP7  
ET,1,45	            ! 趨Ԫ
MP,EX,1,EXX1	    ! ϵĵģ
MP,PRXY,1,NUXY1     ! ϵĲɱ
!****************ģ***************************
*if,kmodel,eq,1,then	! ģ趨ģͲĳ
    L=L	
*else
    L=L/2	! 1/2ģͲ峤ȡ
*endif
CYL4, , ,Di/2, ,Di/2+t, ,L      ! ԲͲ
lsel,s,radius,,Di/2,Di/2+t      ! ѡԲ
LESIZE,all, , ,nhoop,, , , ,1 	! 趨Բʷ
alls                            ! ȫѡ
lsel,u,radius,,Di/2,Di/2+t      ! ȥԲ
LESIZE,all, , ,naxial,, , , ,1 	! 趨ʷ
alls                            ! ȫѡ
wprot,,90                       ! תƽ
vsbw,all                        ! ùƽи
lsel,u,radius,,Di/2,Di/2+t      ! ȥԲ  
lsel,u,length,,L                ! ȥ
LESIZE,all, , ,nt,, , , ,1 	! 趨ȷʷ

vsweep,all	                ! ɨʷ

*if,kmodel,eq,1,then	! ģ趨Լ
    nsel,s,loc,z,0	! ȫģͣԼĳһλ
    nsel,r,loc,x,Di/2+t	    
    D,all,uz,	        ! ȫģͣԼĳһλ
    nsel,s,loc,z,0	! ѡһڵ
    nsel,a,loc,z,L	! ѡһڵ
*else
    asel,s,loc,z,0	! ڰģͣڶԳʩӶԳԼ
    Da,all,symm
    nsel,s,loc,z,L	! ѡһڵ
*endif
csys,1	! 
NROTAT,all	! תѡڵĽڵϵ
D,all,UY	! ʩӻλԼѱ֤ý汣ΪԲ
fini
!**************** ȡ***************************
/SOLU
ANTYPE,STATIC        	    ! ̬
PSTRES,ON	            ! ԤӦӰ
lsel,s,radius,,Di/2+t       ! ѡ뾶ΪԲͲ뾶
asll                        ! ѡ
asel,u,loc,z,0              ! ȥ
asel,u,loc,z,L
asel,u,loc,y,0
cm,aouter,area
SFA,all,1,PRES,p 	    ! Բǵʩѹ
Alls	                    ! ȫѡ
SOLVE	                    ! 
FINI	                    ! ˳
!**************** ȡ***************************
/SOLU	                    ! ½
ANTYPE,BUCKLE 	! ָΪ
BUCOPT,LANB,1 	! Block Lanczosн㣬ȡһģ̬
MXPAND,1                  	! 趨ģ̬չΪ1
SOLVE	! 
FINISH

!*******************************************
/POST1
SET,LAST                       !ȡһغӲ
UPCOORD,0.1,ON                 !趨ʼȱ
SAVE,PERT,DB                   !ʼȱ
FINISH	
RESU,PERT,DB                   !ָʼȱ
!*******************************************
/SOLU
outres,,1                      !ÿһغӲ
ANTYPE,STATIC                  !
NLGEOM,ON                      !η
KBC,0                          !Լ
alls
cmsel,s,aouter
SFA,all,1,PRES,p*2.5           !ֵʩغ
NSUBST,50                      !趨غӲΪ50
ARCLEN,ON,5                    !û
SOLVE
FINI
!*******************************************
/POST1
SET,list,2                     !бȷٽغ

!*******************************************
/POST1
SET,FIRST
PLDISP,1
